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Abstract. With the aim of determining the statistical properties of relativistic tur- 



Or 



bulence and unveiling novel and non-classical features, we present the results of di- 
rect numerical simulations of driven turbulence in an ultrarelativistic hot plasma using 
high-order numerical schemes. We study the statistical properties of flows with average 
Mach number ranging from ~ 0.4 to ~ 1.7 and with average Lorentz factors up to ~ 1.7. 
We find that flow quantities, such as the energy density or the local Lorentz factor, show 
large spatial variance even in the subsonic case as compressibility is enhanced by rela- 
tivistic effects. The velocity field is highly intermittent, but its power-spectrum is found 
to be in good agreement with the predictions of the classical theory of Kolmogorov. 



1. Introduction 

Turbulence is an ubiquitous phenomenon in nature as it plays a fundamental role in 
shaping the dynamics of systems ranging from the mixture of air and oil in a car en- 
gine, up to the rarefied hot plasma composing the intergalactic medium. Relativistic 
hydrodynamics is a fundamental ingredient in the modeling of a number of systems 
characterized by high Lorentz-factor flows, strong gravity or relativistic temperatures. 
Examples include the early Universe, relativistic jets, gam ma-ray-bur sts (GRBs), rela- 
tivistic heavy-ion collisions and core-collapse supernovae (lFontll20(j8 ). 



Despite the importance of relativistic hydrodynamics and the reasonable expecta- 
tion that turbulence is likely to play an important role in many of the systems mentioned 
above, extremely little is known about turbulence in a relativistic regime. For this rea- 
son, the study of relativistic turbulence may be of fundamental importance to develop a 
quantitative description of many astrophysical systems. To this aim, we have performed 
a series of high-order direct numerical simulations of driven relativistic turbulence of a 
hot plasma. 



2. Model and method 

We consider an idealized model of an ultrarelativistic fluid with four-velocity = 
W(l,v'), where W = (1 - v;v') ' 2 is the Lorentz factor and V is the three-velocity in 
units where c - 1. The fluid is modeled as perfect and described by the stress-energy 
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Figure 1 . Left panel: average Lorentz factor as a function of time for the different 
models considered. Note that a quasi-stationary state is reached before t ~ 10 for all 
values of the driving force. Right panel: logarithm of the Lorentz factor on the (y, z) 
plane at the final time of model D. Note the large spatial variations of the Lorentz 
factor with front-like structures. The time-averaged PDFs are shown in the lower left 
corner for the different models considered. 



tensor 

T MV = (P + p)U fl Uy+ Pg/jy, (1) 

where p is the (local-rest-frame) energy density, p is the pressure, u M the four-velocity, 
and gfj V is the spacetime metric, which we take to be the Minkowski one. We evolve 
the equations describing conservation of energy and momentum in the presence of an 
externally imposed Minkowskian force F^, i.e. 

VyT ny _ F (2 ) 

where the forcing term is written as F M = F(0, /'). More specifically, the spatial part of 
the force, f, is a zero-average, solenoidal, random, vector field with a spectral distribu- 
tion which has compact support in the low wavenumber part of the Fourier spectrum. 
Moreover, f, is kept fixed during the evolution and it is th e same for all the models , 



while F is either a constant or a simple function of time (see lRadice & Rezzollal (120121) 
for details). 

The set of relativistic-hydrodynamic equations is closed by the equation of state 
(EOS) p - thus modelling a hot, optically-thick, radiation-pressure dominated 
plasma, such as the electron-positron plasma in a GRB fireball or the matter in the 
radiation-dominated era of the early Universe. The EOS used can be thought as the rel- 
ativistic equivalent of the classical isothermal EOS in that the sound speed is a constant, 
i.e. c\ = 1/3. At the same time, an ultrarelativistic fluid is fundamentally different from 
a classical isothermal fluid. For instance, its "inertia" is entirely determined by the tem- 
perature and the notion of rest-mass density is lost since the latter is minute (or zero 
for a pure photon gas) when compared with the internal one. For these reasons, there is 
no direct classical counterpart of an ultrarelativistic fluid and a relativistic description 
is needed even for small velocities. 

We solve the equations of relativistic hydrodynamics in a 3D periodic domain us - 
ing the high-resolution shock capturing scheme descr ibed i n ( Radice & Rezzollall2012 ). 



In particular, ours is a flux-vector-splitting scheme (lTorolll999h . using the fift h-order 



MP5 reconstruction ( Suresh & Huynh|[l997 ). in local characteristic variables ( Hawke 



200 lj), with a linearized fl ux-split algorithm with entropy and carbuncle fix 



(Radice & Rezzolla 2012). 



Universality and intermittency in relativistic turbulent flows of a hot gas 



3 



3. Basic flow properties 

Our analysis is based on the study of four different models, which we label as A, B, C 
and D, and which differ for the initial amplitude of the driving factor F = 1,2,5 for 
models A-C, and F(f) = 10 + \t for the extreme model D. Each model was evolved 
using three different uniform resolutions of 128 3 , 256 3 and 512 3 grid-zones over the 
same unit lengthscale. As a result, model A is subsonic, model B is transonic and 
models C and D are instead supersonic. The spatial and time-averaged relativistic Mach 
numbers (vW)/(c s W s ) are 0.362, 0.543, 1.003 and 1.759 for our models A, B, C and D, 
while the average Lorentz factors are 1.038, 1.085, 1.278 and 1.732 respectively 

The initial conditions are simple: a constant energy density and a zero-velocity 
field. The forcing term, which is enabled at time t - 0, quickly accelerates the fluid, 
which becomes turbulent. By the time when we start to sample the data, i.e. at ? = 10 
(light-)crossing times, turbulence is fully developed and the flow has reached a station- 
ary state. The evolution is then carried out up to time t = 40, thus providing data for 
15, equally-spaced timeslices over 30 crossing times. As a representative indicator of 
the dynamics of the system, we show in the left panel of Fig. Q] the time evolution of 
the average Lorentz factor for the different models considered. Note that the Lorentz 
factor grows very rapidly during the first few crossing times and then settles to a quasi- 
stationary evolution. Furthermore, the average grows nonlinearly with the increase of 
the driving term, going from (W) ^ 1.04 for the subsonic model A, up to (W) ^ 1.73 
for the most supersonic model D. 

The probability distribution functions (PDFs) of the Lorentz factor are shown in 
the right panel of Fig. [T]for the different models. Clearly, as the forcing is increased, 
the distribution widens, reaching Lorentz factors as large as W 40 (i.e. to speeds 
v - 0.9997). Even in the most "classical" case A, the flow shows patches of fluid 
moving at ultrarelativistic speeds. Also shown in Fig. [T]is the logarithm of the Lorentz 
factor on the (y,z) plane and at t = 40 for model D, highlighting the large spatial 
variations of W and the formation of front-like structures. 



4. Universality 

As customary in studies of turbulence, we have analyzed the power spectrum of the 
velocity field 

E v {k) = ]- f \m\ 2 dk, (3) 
1 J\k\=k 



where A: is a wavenumber three-vector and 

v(k)= \ v(x)e- lnik x dx , (4) 

with V being the three- volume of our computational domain. A number of recent stud- 
ies have analyzed the scaling of the velocity power spectrum in the inertial range, that 
is, in the range in wavenumbers between the leng thscale of the prob lem and the scale 
at which dissipation dominates. More specifically. Ilnoue et all ( 201 1 ) has reported evi- 



dences of a Kolmogorov k 5/3 scaling in a freely-decaying MHD turbulence, but has not 
provided a systematic convergence study of the spectrum. Evidences for a k~ 5 ' 3 scal- 



ing were also found bv lZhangetaH(l2009h . in the case of the kinetic-energy spectrum, 
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Figure 2. Power spectra of the velocity field. Different lines refer to the three 
resolutions used and to the different values of the driving force. The spectra are 
scaled assuming a k~ 5 ^ law. 



which coincides with the ve locity power-spectrum in the incompressible case. Finally, 
Zrake & MacFadyenl ( 2012h has performed a significantly more systematic study for 



driven, transonic, MHD turbulence, but obtained only a very small (if any) coverage of 
the inertial range. 

The time-averaged velocity power spectra computed from our simulations are 
shown in Fig. [2] Different lines refer to the three different resolutions used, 128 3 (dash- 
dotted), 256 3 (dashed) and 512 3 (solid lines), and to the different values of the driving 
force. To highlight the presence and extension of the inertial range, the spectra are 
scaled assuming a kT 5 ^ law, with curves at different resolutions shifted of a factor two 
or four, and nicely overlapping with the high-resolution one in the dissipation region. 
Overall, Fig. |2] convincingly demonstrates the good statistical convergence of our code 
and gi ves a strong support to the idea that the key prediction of the Kolmogorov model 
(K41) (lKolmogorovlll99ll) carries over to the relativistic case. Indeed, not only does the 
velocity spectrum for our subsonic model A shows a region, of about a decade in length, 
where the k~ 5/3 scaling holds, but this continues to be the case even as we increase the 
forcing and enter the regime of relativistic supersonic turbulence with model D. In this 
transition, the velocity spectrum in the inertial range, the range of lengthscales where 
the flow is scale-invariant, is simply "shifted upwards" in a self-similar way, with a 
progressive flattening of the bottleneck region, the bump in the spectrum due to the 
non-linear dissipation introduced by our numerical scheme. Steeper scalings, such as 
the Burger one, k~ 2 , are also clearly incompatible with our data. 

All in all, this is one of our main results: the velocity power spectrum in the inertial 
range is universal, that is, insensitive to relativistic effects, at least in the subsonic and 
mildly supersonic cases. Note that this does not mean that relativistic effects are absent 
or can be neglected when modelling relativistic turbulent flows. 
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Figure 3. Structure function exponents as computed using the ESS technique for 
different models. Also shown in the Figure are the analyti c predictions from two 
classical intermitt ency models: th e She and Leveque (SL) (IShe & Leveque|[T994l) 
and the Boldyrev (iBoldvr ev 2002) models. 



5. Intermittency 

Not all of the information about relativistic turbulent flows is contained in the velocity 
power spectrum. Particularly important in a relativistic context is the intermittency of 
the velocity field, that is, the local appearance of anomalous, short-lived flow features, 
which we have studied by looking at the parallel-structure functions of order p 



SUr) = (\6 r v\ p ), 



8 r v - [v(x + r) - v(x)] ■ - 
r 



(5) 



where r is a vector of length r and the average is over space and time. 

The scaling exponents of the parallel structure functions, i.e. g p s.t. S p (r) ~ r^>, 
have been comput ed up to p - 10 using the extended-self-similarity (ESS) technique 
dBenzi et al.ll 19931) and are summarized in Figure [3] The errors are estimated by com- 
puting the exponents without the ESS or using only the data at the final time. We also 
show the values as computed usi ng the classical K41 t heory, as well as using the esti- 
mates by She and Leveque (SL) dShe & Leveque|[r994l) for incompressible turbulence, 
i.e. 4 = f + 2 - 2(f)^ /3 , and those by Boldyrev dBoldvrevll2002h for Kolmogorov- 
Burgers supersonic turbulence, i.e. t}l = § + 1 - (\) p ^ ■ 
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Not surprisingly, as the flow becomes supersonic, the high-order exponents tend 
to flatten out and be compatible with the Boldyrev scaling, as the most singular veloc- 
ity structures become two-dimensional shock waves, gl, instead, is compatible with 
the She-Leveque model even in the supersonic case. This is consistent with the ob- 
served scaling of the velocity power spectrum, which presents only small intermit- 
tenc y corrections to th e k~ 5 ^ 3 scaling. Previous cl assical studies of weakly compress- 
ible ( Benzi etlul2008h and supersonic turbulence ( Porter et~ail2002h found the scaling 
exponents to be in very good agreement with the ones of the incompressible case and 
to be well described by the SL model. This is very different from what we observe 
even in our subsonic model A, in which the exponents are significantly flatter than in 
the SL model, suggesting a stronger intermittency correction. This deviation is another 
important result of our simulations. 



6. Conclusions 

Using a series of high-order direct numerical simulations of driven relativistic turbu- 
lence in a hot plasma, we have explored the statistical properties of relativistic turbulent 
flows with average Mach numbers ranging from 0.4 to 1 .7 and average Lorentz factors 
up to 1.7. We have found that relativistic effects enhance significantly the intermittency 
of the flow and affect the high-order statistics of the velocity field. Nevertheless, the 
low-order statistics appear to be universal, i.e. independent from the Lorentz factor, and 
in good agreement with the classical Kolmogorov theory. 

Acknowledgments. We thank M.A. Aloy, P. Cerda-Duran, and M. Obergaulinger 
for discussions. The calculations were performed on the clusters at the AEI and on the 
SuperMUC cluster at the LRZ. Partial support comes from the DFG grant SFB/Trans- 
regio 7 and by "CompStar", a Research Networking Programme of the ESF. 



References 

Benzi, R., Biferale, L., Fisher, R., Kadanoff, L., Lamb, D., & Toschi, F. 2008, Phys. Rev. Lett., 
100, 1 

Benzi, R., Ciliberto, S., Tripiccione, R., Baudet, C, Massaioli, F, & Succi, S. 1993, Phys. Rev. 
E, 48, 29 

Boldyrev, S. 2002, The Astrophysical Jour nal, 569, 841 
Font, J. A. 2008, Living Rev. Relativ., 6, 4. 10704 . 2 60871 
Hawke, I. 2001, Ph.D. thesis, University of Cambridge 
Inoue, T, Asano, K., & Ioka, K. 201 1, Astrophys. J., 734, 77 

Kolmogorov, A. N. 1991, Proceedings of the Royal Society A: Mathematical, Physical and 

Engineering Sciences, 434, 9 

Porter, D., Pouquet, A., & Woodward, P. 2002, Phys. Rev. E, 66, 1 

Radice, D., & Rezzolla, L. 2012, Astronomy and Astrophysics. In press. Fl2 6.65 0"2l 

Radice, D., & Rezzolla, L. 2012, Submitted. 1209.29361 

She, Z., & Leveque, E. 1994, Phys. Rev. Lett., 72, 336 

Suresh, A., & Huynh, H. T. 1997, Journal of Computational Physics, 136, 83 

Toro, E. F. 1999, Riemann Solvers and Numerical Methods for Fluid Dynamics (Springer- 

Verlag) 

Zhang, W., MacFadyen, A., & Wang, P. 2009, Astrophys. J. , 692, L40.I0811 . 36381 
Zrake, J., & MacFadyen, A. I. 2012, Astrophys. J., 744, 32. 11108 . l"99Tl 



